NMR and /xSR study of spin correlations in SrZnVO(P04)2 a S=l/2 frustrated 

magnet on a square lattice 
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31 P nuclear and muon spin-lattice relaxation rate measurements in SrZnVO(P04)2, a S = 1/2 
frustrated magnet on a square lattice, are presented. The temperature (T) dependence of the in- 
plane correlation length £ is derived and it is shown that the overall behaviour is analogous to the 
one found for non-frustrated systems but with a reduced spin stiffness. The temperature dependence 
of £ in SrZnVO(P04)2 is compared to the one of other frustrated magnets on a square lattice with 
competing nearest neighbour (Ji) and next-nearest neighbour (J2) exchange couplings and it is 
shown that £ progressively decreases as the ratio J2 / J\ approaches the critical value leading to the 
suppression of long-range magnetic order. In spite of the differences in the functional form of £(T) 
found in different vanadates, it is pointed out that the characteristic energy scale describing spin 
correlations in all those compounds appears to scale as \2J2 + Ji\- 

PACS numbers: 76.60.Es, 76.75.-K, 75.10.Jm, 75.40.Gb 
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I. INTRODUCTION 



The study of quantum magnetism has received a re- 
newed attention after the discovery of high temperature 
superconductivity in the cuprates. In fact, these ma- 
terials have allowed to investigate at the experimental 
level the phase diagram of S=l/2 Heisenberg antiferro- 
magnets on a square lattice with great accuracy^. The 
behaviour of the correlation length has been derived by 
means of neutron scattering experiments 2 - and nuclear 
spin-lattice relaxation rate measurements 3 -^, the form of 
the dynamical susceptibility and, accordingly, the value 
of the scaling exponents which characterize those systems 
have been obtained^ More recently much attention has 
been addressed to the investigation of frustrated square 
lattice (FSL) systems where the frustration is induced 
by a next nearest neighbour (n.n.n.) exchange coupling 
Ji competing with the nearest neighbour (n.n.) one ( J\) 
along the side of the square^. Frustration is expected 
to further enhance quantum fluctuations and to lead to 
the suppression of long-range magnetic order for certain 
values of the ratio r = Jij 3\. In particular, when both 
exchange couplings are antiferromagnetic for r ~ 0.5 a 
spin-liquid ground-state is expected^, while when J\ is 
ferromagnetic (i.e. J± < 0) for r ~ —0.5 a nematic order 
for the two-spin correlation function is envisaged^ The 
Ji — J2 model on a square lattice has received renewed at- 
tention in the last two years when it was realized that the 
parent compounds of the recently discovered iron-based 
superconductors are characterized by comparable n.n. 
and n.n.n. hopping integrals, which may yield competing 
exchange couplings within the square lattice formed by 
iron atoms^. In fact, those materials would represent an 
extension of the J\ — Ji model on a square lattice to itiner- 
ant electron systems. As regards the insulating systems, 
a number of compounds have been recently identified to 
be prototypes of FSL systems and investigated through 



different experimental approaches .^2^— Attempts have 
been made to theoretically understand their high field 
propertie o 10 i 16 and the exchange mechanisms^ Despite 
such a theoretical and experimental progress, a number 
of key questions still have to be addressed. For example, 
the parts of the phase diagram where long-range order 
should be absent have not been studied so far, moreover 
it is not clear how the temperature dependence of the 
in-plane correlation length £ changes with r. Unfortu- 
nately, it is not possible to address this latter aspect by 
means of inelastic neutron scattering experiments since 
only small crystals are available for the prototypes of the 
Ji — J2 model on a square lattice^. Hence, it would be 
worthwhile to find other experimental techniques which 
could allow to determine the effect of frustration on £. 

Here we present an experimental study of the temper- 
ature (T) dependence of the in-plane correlation length 
£, derived by means of nuclear and muon spin-lattice re- 
laxation rates, in SrZnVO(P04)2 (FigHJ) a prototype of 
frustrated magnet on a square lattice with competing 
ferromagnetic n.n. and antiferromagnetic n.n.n. cou- 
plings. It will be shown that £ diverges exponentially on 
cooling with a reduced spin stiffness, possibly scaling as 
I J\ + 2J2I. A comparison with the results previously ob- 
tained by our group on other systems with r < appears 
to qualitatively support this scaling of the spin stiffness, 
even if an accurate description of £ on approaching the 
transition to the columnar ground-state should take into 
account the spin anisotropy and interlayer couplings. 



II. TECHNICAL ASPECTS AND 
EXPERIMENTAL RESULTS 

The synthesis of SrZnVO(P04)2 polycrystalline sam- 
ple was carried out by using the protocol already reported 
in Refsfl^andlTsl DC magnetization (M) measurements 
were performed in order to estimate the superexchange 
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coupling constants for our sample and to check if they are 
consistent with the ones reported in the literature! 10 ' 15 ' 17 
The T-dependence of the static uniform spin suscepti- 
bility x — M/Hr& with H the magnetic field inten- 
sity, was analyzed by fitting the high-T data to Curie- 
Weiss law and to the high-T series expansion.— It was 
found that, Ji = -7.53 ± 0.7 K and J 2 = 8.63 ± 0.6 
K, values which are quite consistent with the ones previ- 
ously reported in the literature ! 10 ' 15 ' 17 In the following, 
in order to better compare SrZnVO(P04)2 to the other 
systems we shall introduce a characteristic energy scale 
Jc = v J i + J 2 — 11-45 K, which provides the magni- 
tude of the exchange couplings. 



with x but with opposite slopes, indicating an opposite 
sign in the hyperfine coupling components (Fig. [2£b)). 




FIG. 1: Projection of SrZnVO(PC>4)2 structure along the c 
axis evidencing the planes containing V 4+ ions. Vanadium 
ions are in green, phosphorus in purple and oxygen ions in 
red. VOs pyramids and PIO4 tetrahedra are also visible. 
The red dotted square shows the 5 = 1/2 square lattice. 

31 P NMR measurements were carried out by using 
standard radiofrequency (RF) pulse sequences. At low 
field, where the full spectrum could be irradiated, the 
NMR powder spectra were obtained from the Fourier 
transform of half of the echo after a7r/2 — te — tt pulse 
sequence. The NMR powder spectrum Fig. [2a) was 
characterized by an asymmetric line shape, quite simi- 
lar to the one found by Nath et ali^I in the isostructural 
Pb2VO(P04)2 compound. The narrow central compo- 
nent is associated with P2 sites lying in between adja- 
cent vanadium layers, while the broader component to 
the PI site (Fig. [IJ which lies within vanadium layers 
and is characterized by a larger hyperfine coupling. At 
high magnetic fields the line broadening prevented the ir- 
radiation of the whole line and the NMR spectrum had to 
be derived either by recording the intensity of the signal 
upon making discrete frequency steps or upon sweeping 
the magnetic field (Fig. Ufa)). The T-dependence of the 
NMR shift AK for the PI site for H || c and H || ab 
was determined by recording the position of the low fre- 
quency (high field) and of the high frequency (low field) 
shoulders of the NMR spectrum, respectively, as a func- 
tion of T. Both quantities were found to scale linearly 
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FIG. 2: (a) Field swept 31 P NMR spectrum is reported for 
RF irradiation at v = 70 MHz. The contribution from 31 P1 
and 31 P2 sites is evidenced and the parts of 31 PI spectra 
corresponding to an orientation of the grains with H || c or 
H || ab are shown, (b) The 31 P NMR shift of the high (green) 
and low (blue) frequency shoulders of 31 P1 NMR spectra, 
corresponding to H || ab and H || c, respectively, is reported 
as a function of the macroscopic spin susceptibility with the 
temperature as an implicit parameter. 



Nuclear spin-lattice relaxation rate 1/T\ was derived 
from the recovery of the nuclear magnetization after a 
saturating pulse sequence. In view of the anisotropy of 
the hyperfine coupling tensor 1/Tj depends on the por- 
tion of the spectrum being irradiated. Since the low- 
frequency (high field) shoulder of the 31 P1 spectrum was 
more separated from the rest of the spectra we have de- 
cided to irradiate just that part of 31 P NMR powder spec- 
trum, corresponding to the crystallites with H || c. From 
now on we shall refer to T\ only for the 31 P1 site and for 
that orientation. The corresponding recovery laws for the 
nuclear magnetization could be nicely fit by a single ex- 
ponential. In Fig. [3]the temperature dependence of 1/Ti 



3 



in the 1.6 K - 100 K range is shown. At high temperature 
1/T\ is flat, then it smoothly decreases and eventually it 
shows a well defined peak at Tc — 2.65 ± 0.02 K, corre- 
sponding to the columnar ordering temperature. Below 
Tc a rapid decrease of 1/Ti is observed. This behaviour 
is very similar to the one reported by Nath et alj^i for 
Pb2VO(P04)2- No significant change in 1/Ti was no- 
ticed upon increasing the magnetic field intensity from 7 
to 35 kGauss (Figj3]) at high T. On the other hand, a 
tiny change has to be expected for T — > T c due to the 
variation of the transition temperature with the field. 21 



characterized by a stretched exponential function A{t) = 
A(0)exp(—(Xt) l3 )r^ with /3 progressively decreasing from 
0.7 to 0.5 upon decreasing the temperature from 30 K 
to Tc- The stretched exponential character of the relax- 
ation can be associated either with a distribution of muon 
sites or with an anisotropic hyperfine coupling, yielding 
to a distribution of relaxation rates in a powder sample. 
Below Tc clear oscillations are observed in zero-field^ 
showing that there is a spontaneous sublattice magneti- 
zation causing a non-zero magnetic field at the muon 
site22. Accordingly the decay of the muon asymmetry 
followed the behaviour typically found for powder sam- 
ples in zero-field 2 ^ 
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FIG. 3: Temperature dependence of 31 P1 nuclear spin-lattice 
relaxation rate in SrZnVO(P04)2. 
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FIG. 4: Temperature dependence of the zero-field muon re- 
laxation rate in SrZnVO(P04)2- 

/iSR measurements were performed at ISIS pulsed 
muon source on MUSR beam line. In zero-field (ZF), 
for T > Tc, the decay of the muon asymmetry was 



A(t) =A 1 e-' 7t cos(j fi B fl t- 



A 2 e 



-At 



B 



(1) 



where 7 M is the muon gyromagnetic ratio, a the decay 
rate of the oscillating part, mostly due to a inhomoge- 
neous distribution of the local field at the muon sites, 
while B is a constant background arising from the sam- 
ple environment. 




FIG. 5: The local field at the muon, normalized to its low tem- 
perature value Bfj, — 186 ± 3 Gauss, is reported as a function 
of T/T c , with T c = 2.65 K. The low-temperature blue dotted 
line shows the behaviour expected for a power law reduction 
of B M (T) = B M (0)(1 - aT 2 ). The high temperature red dot- 
ted line shows the critical behaviour for a critical exponent 
P = 0.235, expected for a 2D XY model. 

The temperature dependence of A and of B^ derived 
from the fit of the asymmetry with the aforementioned 
expressions is reported in Fig. 0]and Fig. respectively. 



III. DISCUSSION 

First we shall consider the temperature dependence 
of the order parameter, as derived from zero-field /iSR 
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measurements. The local field at the muon can be written 



J2^<Si>=A» ff \<S> 



(2) 



where is the hyperfine coupling between the muon 
and the i-th V 4+ spin. Since the magnitude of all V 4+ 
spins | < S > | is expected to be the same, the local 
field at the muon can be written in terms of an effective 
total hyperfine coupling A^tt times | < S > |. Hence, 
the T-dependence of the local field at the muon gives 
directly the one of the sublattice magnetization. As re- 
gards, the critical behaviour of the order parameter, here 
we only remark that the critical exponent j3 can be con- 
sistent with the one expected for finite 2DXY systems 
(j3 = 0.235)2^ found in other similar vanadatesj 21 ' 22 i 25 
Nevertheless, the accuracy of the experimental points 
does not allow one to give a definite answer in this re- 
spect. On the other hand, the low temperature behaviour 
of provides information on the dispersion relation for 
the spin waves excitations. The reduction of the low- 
temperature sublattice magnetization is consistent with 
a power-law | < S > |(T) ~ T n , with n = 1.8 ± 0.4. Al- 
though this value would be consistent with the dispersion 
relation for nearly two-dimensional antiferromagnets^ it 
is difficult to give a precise statement in view of the ex- 
perimental uncertainty. Nevertheless, as it will be shown 
at the end of this section, also the T-dependence of 1/Ti 
seems to support the 2D character of the spin wave exci- 
tations, with a quasi-linear magnon dispersion. 

Now we turn to the discussion of the temperature 
dependence of P nuclear spin- lattice relaxation rate, 
which allows to derive information on the low-energy dy- 
namics and on the spin correlations. In the case of a 
magnetic relaxation process driven by electron spin fluc- 
tuations 1/Ti can be written 2 ^ 



(3) 



where 7 is the nuclear gyromagnetic ratio, |^4 q | 2 the form 
factor describing the hyperfine coupling with spin excita- 
tions at wave-vector q and S a ,a(<l, ^l) (a = x,y,z) the 
component of the dynamical structure factor at the Lar- 
mor frequency u^. The _L subscript indicates that one 
should consider the components of the fluctuating hyper- 
fine field perpendicular to the quantization axis, given by 
the direction of the static external field. By using scaling 
arguments it is possible to write the dynamical structure 
factor in terms of the in-plane correlation length £ (in 
lattice units hereafter) and establish a one to one rela- 
tionship between 1/T± and £. This procedure has proven 
to be very useful to study the temperature dependence 
of the correlation length in the cuprates which are proto- 
types of two-dimensional S — 1/2 Heisenberg antifcrro- 
magnets on a square lattice^ and to determine the value 
of the dynamical scaling exponent z = 1. Given the 
similarity between the cuprates and the vanadates under 



investigation, we will use the same approach here to de- 
rive £ from P 1/Ti, assuming 2 = 1. Accordingly, one 
can write^ 



1 2 



s(s+i)_ e+2 m 2 v^ 1 



4tt 2 



bz i + e(q-Qc) 2 



(4) 



where (3 is a normalization factor which allows to pre- 
serve the spin sum rule and Qc is the columnar criti- 
cal wave- vector, uje — {JckB/K)\j2nS(S + l)/3 is the 
Heisenberg exchange frequency^ with n — 4 the number 
of n.n. and of next n.n.. 

In order to establish a one to one relationship between 
1/Tl and £ from the above equation one has first to de- 
rive the hyperfine coupling tensor components and the 
form factor for 31 P1 site. The components of the hy- 
perfine tensor can be derived from the plot of the shift 
AK for the magnetic field along the a direction vs. the 
molar macroscopic susceptibility \ (Fig- 2(b)), which is 
assumed to be isotropic. Then one can write 



AK a = 



9VbNa 



(5) 



with [Ib the Bohr magneton and Na the Avogardro's 
number. 




FIG. 6: The form factors for 31 P1 site are reported as a func- 
tion of the in-plane components (q x and q y ) of the wave- vector 
of the spin excitations for SrZnVO(P04)2- 

For SrZnVO(P0 4 ) 2 one finds A cc = -4300±190 Gauss 
and A aa ~ Abb = 2360 ± 200 Gauss. In fact, since the 
measurements on powders did not allow one to discern 
between AK a and AKb we have assumed A aa ~ Abb- 
In Pb2VO(P04)2, on the other hand, a small difference 
between those two components is observed^ The total 
hyperfine tensor is the sum of a transferred term A* and 
of a dipolar term A d . The latter one can be calculated 
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FIG. 7: 31 P1 nuclear spin-lattice relaxation rate in 
SrZnVO(P04)2 is reported as a function of the in-plane cor- 
relation length (in lattice units) according to Eq(?]in the text. 



on the basis of lattice sums while the former one is as- 
sumed to be the sum of four equal terms arising from the 
hyperfme coupling between 31 P1 nuclei and the four n.n. 
V 4+ spins. Hence the contribution to the transferred hy- 
perfme term is simply given by A* = (A — A d )/4. Now 
that both the transferred and dipolar coupling between 
31 PI nucleus and each V 4+ spin are known, it is possible 
to derive the hyperfme form factor. The form factor of 
SrZnVO(P04)2 is reported in FigJHl It is noticed that, 
owing to the symmetry position of PI site, the form fac- 
tor shows a non-vanishing minimum at Qc — (i 7 ^ 0) and 
(0, ±7r). This explains why 1/Ti progressively decreases 
as the system gets more correlated upon decreasing tem- 
perature (FigI2J) and only when the correlation length is 
sufficiently large 1/Ti increases again. In fact, from EqQ] 
it is possible to derive numerically the behaviour of 1 jT\ 
vs. £ (Fig. [7]) and one finds a minimum for £ ~ 5 lattice 
steps. Accordingly, from the experimental data reported 
in Fig. [3] it is now possible to derive quantitatively the 
temperature dependence of £ in SrZnVO(P04)2- 

In order to derive the temperature dependence of £ 
from A(T) one has first to subtract the T-independent 
dipolar contribution from the raw data in Fig. 4 and then 
proceed in the same way as it was done for 1/Tl. How- 
ever, here the form factor cannot be determined since the 
muon site and hyperfine couplings are unknown. It is no- 
ticed that, at variance with PI site, the muon should not 
be in a position symmetrical with respect to the neigh- 
bouring V 4+ ions, since A continues to diverge upon cool- 
ing. Thus, if no filtering effect due to the form factor is 
present one can safely write A ~ £, for £ 3> 1,1 Namely, 
the temperature dependence of A directly gives the one 
of £ although, unlike 31 P 1/Tl, A does not allow to es- 
timate quantitatively £. Nevertheless, by matching the 
£(T) derived from A(T) with the one quantitatively de- 
rived from 1/Ti over the same T range, it is possible to 
use also A(T) data to estimate quantitatively £(T). 

In Fig|5] we report the temperature dependence of £ 




FIG. 8: The temperature dependence of the in-plane corre- 
lation length (in lattice units) derived from A and 1/Ti data 
is reported as a function of Jc/T. The solid line shows the 
behaviour expected for a spin-stiffness p s = 0.79 ± 0.05 x 
1.15Jc/2tt. 



10 ■ 



r=-1. 
a 



r=-1.15 



r=-0.9 



# B 

a b h 
_ a 



3 



9 



n Pb 2 vo(P0 4 ) 2 

B SrZnVO(P0 4 ) 2 
3 BaCdVO(P0 4 ) 



— 1 1 1 1 1 1 1 1 — 

2 3 4 5 6 

J c /T 

FIG. 9: The temperature dependence in-plane correlation 
length derived from A is reported vs. Jc/T for compared 
to the one derived for SrZnVO(P0 4 ) 2 (Jc = 11-45 K), 
BaCdVO(P0 4 ) 2 (Jc = 4.8 K) 14 and Pb 2 VO(P0 4 ) 2 (Jc = 
10.7 K)^. The dotted arrow points out that upon decreasing 
| J 2 / Ji | the correlation length increases less rapidly on cooling. 



derived by means of 31 P1 1/Ti and the one obtained by 
means of A(T) in a temperature range where £ is suffi- 
ciently large so that either Eq|4] apply or A ~ £. One 
notices an overall good agreement in the behaviour of 
£(T) derived through both methods, moreover it is no- 
ticed that £ diverges exponentially on decreasing tem- 
perature. For a non frustrated S = 1/2 Heisenberg 
antiferromagnet on a square lattice one would expect 
that £(T) ~ exp{2-Kp s /T)/{T + Anps)^ with p s the 
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spin stiffness, which for non frustrated systems turns out 
p s ~ 1.15Jc/27r. Here we find that the behaviour of 
£(T) is the same but with a reduced effective spin stiff- 
ness constant. In fact, the data in FigJ5]can be nicely fit 
with p s = 0.79 ± 0.05 x 1.15 J c /2tt = 1.66 ± 0.1 K. 

Recently Hartel et al^ have calculated the temper- 
ature dependence of £(T) for J 2 < 0.44|Ji| (i.e. r > 
—0.44), namely for the part of the phase diagram ad- 
jacent to the one experimentally investigated here and 
in Refl22l They found that £(T) diverges exponentially 
with decreasing temperature with an effective spin stiff- 
ness p s ~ — {J i + 2J2)/8, which vanishes on approaching 
r ~ —0.5, namely the region with no long-range mag- 
netic order SrZnVO(P04)2, however, is characterized 
by r ~ —1.15 and it is not clear if the previous expres- 
sion for the spin stiffness can still be used. Neverthe- 
less, if one considers that also for the compounds with 
r < —0.5 an analogous expression p s ~ + (Ji + 2J 2 )/8 
could hold, one would derive for SrZnVO(P04)2 an ef- 
fective spin stiffness p s — 1.23 ± 0.16 K, a value which is 
close to the one experimentally determined here (FigJS]). 
Accordingly, in the absence of a theoretical calculation, 
one would be tempted to argue that on both sides of 
the critical point around r ~ —0.5 the correlation length 
diverges exponentially on cooling with an effective spin 
stiffness p s ~ |(Ji + 2J 2 )/8|, progressively vanishing as 
r^0.5. 



£ — s- 1, £ is the same in all compounds. The correspond- 
ing data are reported in Fig. [§1 One notices that indeed 
£ decreases as r — > —0.5, however, it is also noticed that 
while for SrZnVO(P04)2 the correlation length diverges 
exponentially over a wide T range , this is not the case 
for the other two systems. In fact, it has been pointed 
out that the behaviour of Pb2VO(P04)2 is more charac- 
teristic of a 2D XY system, 22 - while in BaCdVO(P0 4 ) 2 
possibly nematic correlations appear, leading to a loga- 
rithmic increase of £ on cooling. 22 Moreover, deviations 
associated with the critical behaviour are observed on 
approaching Tq- Hence, it appears that although in gen- 
eral the system becomes less correlated as r — > —0.5 the 
correct analytical form of the correlation function is not 
simply exponential and that details taking into account 
the presence of a possible XY character or of nematic cor- 
relations should be considered. Nevertheless, it is inter- 
esting to observe that in spite of the functional form, the 
characteristic energy scale describing the growth of the 
in-plane correlation length appears to scale as Ji+2 Ji far 
from Tq- In fact, if one now plots the £ data when £ 3> 1 
for the different compounds as a function of ( J\ +2 J-i) /T, 
one observes a reasonable overlap between the data of the 
different compounds until when the XY character or the 
inter-layer coupling do not give rise to a critical enhance- 
ment of the correlations on approaching Tc- It is noticed 
that also in this plot the magnitude of £ in Pb2VO(P04)2 
and BaCdVO(P04)2 has been rescaled in order to match 
the one quantitatively derived for SrZnVO(P04)2- Hence 
a definite answer on the validity of the scaling would re- 
quire an independent quantitative estimate of £ also for 
those two compounds. 



FIG. 10: The in-plane correlation length derived from A(T) 
is reported as a function of (J a + 2J 2 )/T for SrZnVO(P0 4 ) 2 
((Ji + 2J 2 ) = 9.73 K), BaCdVO(P0 4 ) 2 ((Ji + 2J 2 ) = 2.8 
K)±^i and Pb 2 VO(P0 4 ) 2 ((Ji + 2J 2 ) = 13.7 K)M 

Now, it is rather interesting to compare the behaviour 
of £ in SrZnVO(P0 4 ) 2 with the one in BaCdVO(P0 4 ) 2 
(r ~ -0.9 14 ) and in Pb 2 VO(P0 4 ) 2 (r ^ -1.8 28 ), de- 
rived from the temperature dependence of A(T)j22 Since 
in the latter two compounds it was not possible to de- 
termine the absolute value of £ we assumed that A ~ £ 
and rescaled the values of A so that for T ~ Jc/2, when 
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FIG. 11: The temperature dependence of 31 P1 1/Ti in 
SrZnVO(P04) 2 in the columnar ground-state is reported. 
The solid line shows the best fit according to a power-law 
behaviour l/Ti ~ T b with b = 1.9 ± 0.3. 

Finally, we shall discuss the behaviour of the spin- 
lattice relaxation rate in the magnetically ordered colum- 
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nar phase. Below Tq one observes a marked decrease of 
1/Ti which in the low-temperature limit should be as- 
cribed to the vanishing of the two-magnon Raman re- 
laxation processes^ If the gap in the magnon disper- 
sion curve is negligible one would expect a power-law be- 
haviour of 1/Ti with a power law exponent depending on 
the magnetic lattice dimensionality and on the analytical 
form of the magnon dispersion curved In case of a linear 
dispersion curve, neglecting the presence of a gap in the 
spin-wave dispersion, for a quasi-2D system one would 
expect 1/Ti - T 2 . Here we find that 1/Ti ~ T b with 
b = 1.9±0.3, in reasonable agreement with the theoreti- 
cal expectations and with the behaviour of the sublattice 
magnetization derived from /i + SR measurements. 

IV. CONCLUSIONS 

In conclusion we have determined quantitatively the 
temperature dependence of the in-plane correlation 



length £ in SrZnVO(P0 4 ) 2 , a frustrated 5 = 1/2 mag- 
net on a square lattice with r ~ —1.15, by means of 
nuclear and muon spin-lattice relaxation rate measure- 
ments. It has been shown that £ diverges exponentially 
on cooling with a reduced spin stiffness which appears to 
roughly scale as | J\ + 2 J 2 \ . A comparison with the results 
previously obtained by our group on other systems with 
r < appears to support this scaling of the spin stiffness 
even if an accurate description of £ on approaching the 
transition to the columnar ground-state should take into 
account the spin anisotropy and interlayer couplings. 
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